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The two-dimensional model of itinerant electrons coupled to an anti-ferromagnetic order parame- 
ter is considered. In the mean field solution the Fermi surface undergoes reconstruction, and breaks 
into disconnected "pockets" . We have studied the effect of the thermal fluctuations of the order 
parameter on the spectral density in such system. These fluctuations lead to a finite width of the 
spectral line scaling linearly with temperature. Due to the thermal fluctuations the quasi-particle 
spectral weight is transfered into a magnetic Brillouin zone. This can be interpreted as restoration 
of "arcs" of the non-interacting Fermi surface. 
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I. INTRODUCTION 

Fluctuations play a prominent role in systems of re- 
duced dimensionalities leading to a complete or partial 
suppression of long range order and rendering mean field 
approximation inapplicable. Calculation of correlation 
functions then becomes an arduous task. The question 
is whether in the absence of true long range order these 
functions display features of the ordered state and if yes, 
to what extent. The correlation function we are con- 
cerned with in this paper is the single electron spec- 
tral function. Measurements of this function constitute 
the most powerful experimental probes in the physics of 
strongly correlated systems. We discuss the situation 
when the system is close to being antiferromagnetically 
ordered and study the effect of thermal order parameter 
fluctuations. 

There is a considerable literature addressing the influ- 
ence of quantum fluctuations (we refer the reader to [l| 
and [|[ which also provide references to the related pa- 
pers). However, at finite temperatures one has to take 
into account thermal fluctuations which brings specific 
problems. In our previous publication [H we consid- 
ered the influence of thermal fluctuations on the spectral 
function of a two dimensional (2D) superconductor. 2D 
superconductors have quasi long range order such that 
phase fluctuations are critical in the entire temperature 
region below T c . We have found that at least as far as 
the thermal fluctuations were concerned, the frequently 
used approach based on the conversion of this problem 
into a gauge theory turned out to be inadequate. The 
latter approach includes a gauge transformation of the 
fermion fields with a subsequent attempts to treat the 
problem as a gauge field theory one (as, for instance, 
in [2|). The difficulty comes from the fact that the re- 
sulting calculational scheme contains strong ultraviolet 
divergencies and hence requires knowledge about states 
located far from the Fermi surface. As an alternative we 
have suggested the direct perturbation expansion in the 
order parameter. This procedure contains only infrared 
divergencies and is therefore universal. 



In the present paper we apply the approach of Q to cal- 
culate the spectral function in the presence of a fluctuat- 
ing commensurate Spin Density Wave (SDW). This prob- 
lem has a potential relevance to the problem of cuprates. 
There is a significant experimental evidence in favor of 
Fermi surface reconstruction taking place in the under- 
doped phase of the copper oxide superconductors (see, 
for example, [H],[5[). On the other hand, it is still un- 
clear whether one needs a real long range order to observe 
such reconstruction or a short range one will suffice. In 
the present paper we will address this question in the 
context of the spectral function. We consider only clas- 
sical (thermal) fluctuations of the order parameter. This 
is more than an academic excercise since thermally fluc- 
tuating SDW phase has been suggested to occupy a part 
of the cuprate phase diagram in the strongly underdoped 
regime 



II. DESCRIPTION OF THE MODEL AND THE 
RESULTS. 

We consider a popular spin-fermion model in two- 
dimensions where electrons interact with a commensu- 
rate SDW [l[ . The antiferromagnetic ordering open gaps 
at the points of the Fermi surface (FS) connected by the 
vector of anti-ferromagnetic fluctuation Q = (7r,7r), see 
Fig. [T] As a result the FS undergoes reconstruction into 
disconnected pockets. The Hamiltonian for quasiparti- 
cles located near two FS points connected by Q is 

H = H £( fc )^L^ + J Y1 s^l+Q^^kp , (i) 

ka k 

where k is a momentum vector in the Brillouin zone, and 
a is the spin index. The kinetic energy close to the "hot 
spots" can be approximated as 



(2) 
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i=l,2 ka 



with the index i enumerating two subbands related by 
the vector Q. The sum in Eq. © is over small mo- 
menta and two subbands can be combined together (see 
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inset in Fig. [TJ. We assume that there is no nesting; to 
simplify the calculations we consider the case when the 
corresponding Fermi velocities are perpendicular to each 
other: v\ — vx, and t>2 = vy. At the mean field level the 
spectrum is determined by the equation 



oj 2 — v(k x + k y )uj + v 2 k x k y — J 2 = 
with the solutions 



(3) 



wi,2 = v(k x + k y )/2 ± ^j[v(k x - k y )/2} 2 + J 2 (4) 
signifying tips of electron- and hole-like Fermi pockets. 




FIG. 1: Formation of a Fermi pockets (thick solid line) in 
the mean field Fermi surface reconstraction caused by the 
SDW ordering. The bare Fermi surface (thin solid line) is 
not nested. The dashed line is the magnetic Brillouin zone 
boundary. Two subbands with i = 1,2 are connected by the 
anti-ferromagnetic wave vector Q = (n,n). Fermi velocities 
vi and V2 are assumed to be perpendicular. Inset shows two 
subbands brought together. 

Fluctuations of the SDW order parameter transform 
rigid energy gaps into pseudogaps and smear the sharp 
peaks in the spectral function. Below we will study this 
process in detail. 

It has been demonstrated in [l[ that the feedback from 
the quasiparticles onto the spin Hamiltonian makes sig- 
nificant changes in the spin dynamics, but does not affect 
zero frequency modes. Since we consider only classical 
(thermal) fluctuations, such feedback will be neglected. 

For the sake of simplicity we also assume that the vec- 
tor of spin polarization lies either (i) in the XY-j>\ane or 
(ii) directed along the z axis. In both cases the transition 
occurs at finite temperature. In the first case the order 
parameter has U(l) symmetry and the transition is of the 
Kosterlitz-Thouless type. The spin fluctuations below T c 
are critical. This power law behavior will also hold above 
T c , though only at distances smaller than the correlation 
length £. However, since the latter length is exponentially 
large in (T — T c ), there is a range of temperatures T > T c 
and energies where the obtained expressions for the spec- 
tral function remain valid. In all this region where the 
magnetic correlation length is either infinite (T < T c ) 



or exponentially large, the order parameter fluctuations 
can be considered as classical. This is essential for our 
approach. In the second case the order parameter has 
Z-i symmetry and below T c it acquires finite expectation 
value. Hence T <T C region is trivial and we will be con- 
cerned only with T > T c region. The correlation length 
in this region is £ ~ (T — T c ) ; to neglect quantum fluc- 
tuations we need it to be much larger than T^ 1 meaning 
that we need to stay close to T c . 

We will start with the easy plane anisotropy case; the 
easy axis case can be obtained as a simple generalization. 
The fluctuating order parameter is staggered magnetiza- 
tion S, it lies in a plane and forms an angle <fi with the 
fixed direction in the plane. Under the assumptions de- 
scribed above the quasiparticle Lagrangian is simplified: 

i,a, 

+ e^i^cr^gV^ + c.c. , (5) 

otfl 

where £i(fc) = Vik. We assume that the free energy for 
the classical phase field <fi is Gaussian: 



F 
T 



Ps_ 

2T 



dxdy [{d x 4>) 2 + {d v <t>f 



(6) 



Now the problem looks similar to the one of the thermal 
fluctuations in superconductors considered in our previ- 
ous paper 0]. 

To be definite we consider the propagator of the spin- 
up particles. The spectral weight reaches its maximal 
value in the vicinity of the bare mass shell, ui ~ k x 
and also close to the mass shell of the spin-down particle 
oj ~ k y (the shadow mass shell). These two regions form 
two complementary parts of the Fermi pocket. As the 
spectral weight is small at the magnetic Brillouin zone 
boundary, k x ~ k y we have studied the Green function 
separately at uj ~ k x and uj ~ k y . 

The summary of our results is as follows. At the mass 
shell we got the following expression for the Green func- 
tion, 



! x 2rfq^r 2 (2 - 2d) J* x 

G -°mf + ( 777. 1 \\4-4d G mf m 



(-i(0J-ky))^ 

2dia 6d T 2 (2 - 2d)r(l - 2d) J 6 

(-iicj-ky))^ 



G 



ad 



(7) 



where d = T/ATrp s is the scaling dimension of the order 
parameter, a is the lattice constant, and the mean field 
Green function is 



G mf = w - A* 



iJ 2 a 2 T(l - 2d) 

(-i(u-ky)y-™ 



(8) 



We would like to emphasise that the expression does 
not rely on the smallness of the parameter d < 1/2. 



At the shadow mass shell, w ~ k y we get 



G 



J 2 e l7Td T(l-2d) 

(w - k x ) 2 



U 2 a 2d T{l - 2d) 



(9) 



Our analytical results, are presented graphycally in 
Fig. [2] The area of validity of these expressions is con- 
trolled by the energy scale 



Tk = J (J a) 



d/(i-d) 



(10) 



Equation ([7]) is valid for \k y \ > Tk/v. Equation © is 
valid for \k x \> Tk/v. 

The remaining part of the paper contains a deriva- 
tion of the results, Eqs. d7j) , (|9j) . Following the approach 
of [|[ we develop a perturbation theory in the coupling 
constant J. Though this perturbation theory is free of 
ultra-violet singularities it contains infra-red singularities 
at w ~ k x (y) which we sum up. 



III. BEHAVIOR AT THE MASS SHELL: uj ~ k x 

In this section we derive the result ([7]). It is useful to 
consider the self-energy E u (fc) defined in the standard 
way by the Dyson equation 



G u {k) = [u - k x - £„(*:)]" 



(11) 



As is shown below, the self-energy is a regular function 
of frequency at the mass shell, and only weakly logarith- 
mically non-analityc in the coupling constant. 
The J 2 contribution to the self-energy is 



4 2) (*0 = 



-i( Ja d yT{l - 2d) 

(-i( w - k y )y- 



2d 



(12) 



We notice that in the limit of infinite phase stiffness, d — 
0, Eq. (TTZj) reproduces the mean field spectrum, Eq. Q, 
as expected. Note that the self energy Eq. (|12|) is regular 
at the mass shell. This, however, is not the case for higher 
order contribution. In fact, the fourth order contribution 
has a weak logarithmic singularity (see Add. IA lj) . 



„ ,-i(Ja d ) 4 r 2 (2-2rf) , 



where 



iO 

T(T 



(13) 



(14) 



The aforementioned analyticity of the self-energy at the 
mass shell is restored once the leading on-shell singulari- 
ties in all orders in J 2 are summed up. The reminder of 
the present section is devoted to this task. 



A w=0 (k) 
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FIG. 2: (color online) The spectral density at the Fermi sur- 
face, A„ = o(fc), as given by Eqs. (0 and ((9]) at k y > k x and 
k x > k y respectively. These two graphs are separated by the 
area where the presented derivation ceases to be valid. The 
parameter d is (a) d = 0.08, and (b) d = 0.2. 



Using the expressions for the bare (retarded) Green 
functions 



iG i(2)( w > r ) = e { r x(y))5{r y (x))e L 



(15) 



we write the self-energy at the order J 2n with n > 3 in 
the form (see Fig. 3) 

poo POO 

S^ 2n) (fc) =i(-iJ) 2 "7 dx n e i{uJ - k * )x ™ d^e^"^^" 



(16) 



11/ dxtYl dy. 



x C (2n) (ri,...,r„; Pl ,...,p n ). 
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where n = (a^yj-i), p t = (xi,yi), x x = y = (see 
Fig. |3]) and the cumulant in the last line 

C (2n) (ri, . . .,r n ;pi, . . . ,p n ) = 5^ n 
n-l 

- ^ *l,i*i+l,n + • • • + (-l)"<5i,i • ■ ■ (17) 

1=1 

is expressed in terms of averages of the exponents of the 
phase fields: 



H,i+l 



- /„10(T- i ) + ... + 10(l-i+l)p-10(Pi)---10(Pi + i) 



(18) 
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FIG. 3: (color online) Graphical representation of the self- 
energy correction of the order 2n in coupling constant, 
Eq. (|16p . Solid vertical and horizontal lines represent seg- 
ments of a real space electron trajectory for (a) particle close 
to the mass shell, tj « k x , (b) particle close to the shadow 
mass shell, cj w Incoming (blue) and outgoing (red) ar- 
rowed skew dashed lines represent exponential factors e 1 ^^^ 
and e~ 1 't'( p 3' respectively. 

The latter average with free energy, Eq. © is well 
known, 



S hn =a 2<M»-2), 



~ /i ~ ft 

{I ., i '•/ II, ; i |Pi-pjl 



n 



■Pi I 



2<i 



(19) 



where tilde over the product sign excludes i = j. 

In what follows we sum up the most singular terms 
in expansion (|16p . Our solution is based on the physical 
idea that at the mass shell the horizontal segments in the 
staircase diagram in Fig. ISl^a) are parametrically longer 
then the vertical ones. In other words, the mass shell 
singularity comes from the integration region yi <C Xj. 
Accordingly, we introduce new variables £j, 



expansion gives 

n 

1=1 

x ' [xl + ( yn ^-y r} d K + (yn-yi) 2 } d 

_K + (Vn - yo) 2 } d [xl + (Vn-l ~ Vl) 2 } d 

n 

^-2da 2dn (x n n^ n Y[^ 2d . ( 21 ) 

i=l 

We substitute Eqs. fH) and dHJ in Eq. (JTBJ), and per- 
form integrations over x n . We now analyze the remaining 
integrations over £jS, 



E< 2 ") =2d 
where 



( _ i)2 » + l (Jarf) 2n r[n(2 _ 2rf) _ 1] 

(n - 2)!(-i(w - fc y ))"( 2 - M )-! M ' [ ' 



w=/ n^ ,. \ , ; 2 . r^-i - ( 23 ) 



o ti +& + ...+&.) 



We notice that for n = 2,3 the integral in Eq. ((23|) di- 
verges at the upper limit. In this case the expansion in 
Eq. (f2"Tj) is not justified. These values of n have to be 
treated separately, (see App. |A1 for details). For n = 2 
the most singular part is given in Eq. (|13|) and for n = 3 
we obtain, (see App. I A 2| ). 



For n > 3 the integral in Eq. (|2"3"|) is 

J„(a) - (1 -2d) 2 a 



, 3 _„r"(l-2d)r(n-3) 



I>(2 - 2d) - 1] 
The contributions of the orders n > 3 give 

^ S ( 2 «).2^-- i (^) 4r2 ( 2 - M ) 

n>3 



(25) 



2d- 



(-i( W -fc,))3-4d 

-i(Ja d ) 6 r(l - 2d)r 2 (2 - 2d) 



log(l + x) 



(-i( w -fc v ))»-<W 



[l-log(l + z)], (26) 



where 



q- 1 (Ja d ) 2 r(l-2d) 
(-i(w - fc v ))( 2 - a «») 



(27) 



The sum of contributions (fTS"]). {T3J), (JHD and ([211 yields 
the final result Eq. ©. 



2M - 2/n-l = £,nXn, ■ ■ ■ , Vl ~ 2/0 = Cl^n (20) 

and argue that the important domain of integration is 
<C 1. We expand the cumulant (fT7|) in powers of £jS 
and retain the lowest power term to get the most singular 
contribution. It can be shown by inspection that this 



IV. GREEN FUNCTION AT THE SHADOW 
MASS SHELL, u ~ k y . 

In this section we turn to the analysis of the behavior 
of the Green function at the "shadow side" of the pocket, 
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u> ~ k y . As this region is separated from the mass shell, 

instead of the self-energy it is more convenient to study 

the Green function itself. It is also convenient to discuss 

— r /q\ "I — 2 

the amputated propagator, £ w (fe) = G u] {k) G u (k) 

To the second order in J, £ u (fc) is given by Eq. (|T2"j). 
and is singular at tu = k y . At the next, fourth order, the 
correction (see App. [B| has stronger singularity 



g(4) 



i( Ja d ) 4 r(l - 2d)T(2 - 2d) 2 _ 2d 
(-i(w - M) 3 " 4d " 



(28) 



In what follows we resum the most singular terms in 
the expansion of £ to all orders in J. For n > 2 we 
introduce new variables, 



■x 2 -xi,... 

X - X n - 2 , 



(29) 



and integrating over y^s variables we write the correction 
of order (2n) to the amputated Green function as 



{2n) (-y-y«") 



(-i(w - fc^)) 2 '^ 1 -' 1 )- 1 
r[2n(l-d)-l] r 

X f 7T1 In{a). 

(n-1)! 



(30) 



where the remaining integrals 



4(a) = n/ * 



-2d 



t-2d 
' Sn-1 



r^pi (3i) 



(a-i+a + .-.+Cn-i) 2 ^ 1 -^ 



are convergent at the upper limit, and can be evaluated 
as 



I n (a) = a 



n-2d 



T n - 1 {l-2d)T{n-2d) 



r(2n(l - d) - 1) 
As a result we obtain for the singular part 

i{-Ua d ) 2n a n - 2d T n - 1 (l-2d)T{n-2d) 



h(2n) _ 



(n - l)!(-i(w - k x )) 2 < 1 - d )- 



(32) 



(33) 



For n — 1 the last expression reduces to the second order 
corrections, Eq. (|12j) . The sum of singular contributions 
Eq. (33} yields the result of Eq. ©. 



V. CONCLUSIONS 

First we would like to make a remark about the easy 
axis anisotropy regime where the phase transition is in 
the Ising model universality class. From our calculations 
it is easy to see that as far as the singular terms are 
concerned, the results remain unchanged provided one 
considers only one particular value for the scaling dimen- 
sion: d = 1/8. This is despite the fact that multipoint 
correlation functions of the Ising model order parame- 
ter fields are more complicated than ([TO)) . However, the 



singularities are determined by more simple correlators, 
namely by the diagrams where pairs of the operators are 
very close to each other (see Fig. |3J resulting in a fusion 
of two order parameter fields. In the Ising model such 
fusion generates the energy density operator and in the 
XY-model it generates the gradient of <fi field. Both op- 
erators have the same multi-point correlation functions. 

Now we can discuss the results. They are well illus- 
trated by Figs. Ill II The region of applicability of our 
calculations involves the energy scale Tk — J(Ja) d ^ d ~ 1 \ 
The result at the mass shell is valid for \k y \ > Tr- and the 
result at the shadow mass shell holds at \k x \ > Tk- With 
increasing temperature the spectral weight is transferred 
towards the bare Fermi surface and the shadow band fea- 
ture quickly fades away as is clearly seen on figure ^ 
where the darker areas correspond to large values of the 
spectral density. 

Although our model does not include all the features 
ascribed to the cuprates, the results obtained may serve 
as a good qualitative guide to the problem. For instance, 
we see that critical thermal fluctuations give rise to the 
characteristic linear temperature dependence ~ T of the 
spectral peak width. This is an indication that such fluc- 
tuations are responsible of this feature in the cuprates. 
Our results demonstrate that with the rise of temperature 
the renormalized mass shell identified as the maximum 
intensity line in Fig. U approaches the bare Fermi surface, 
while the peak becomes rather incoherent. The back- 
side of the Fermi pockets fades away so that the pockets 
now look like arcs. These effects are qualitatively simi- 
lar to that of the quantum fluctuations studied in Ref. @ 
though the intensity of quantum fluctuations is regulated 
not by temperature, but by the interactions. 

In summary, we have presented systematic study of 
the thermal fluctuations effects in two-dimensional sys- 
tem of electrons in interaction with SDW order param- 
eter. In particular the spectral density has been found 
to be strongly sensitive to these fluctuations. Fluctua- 
tions tend to restore the non-interacting Fermi surface 
topology thus overriding the effects of the SDW order. 
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Appendix A: Perturbation theory at the mass shell. 

In the present appendix we evaluate the most singular 
corrections to the self-energy at the mass shell, us ~ k x . 



G 




FIG. 4: False color plot representing the spectral density at 
the Fermi surface, A u= o(k), as given by Eqs. ([7| and at 
k y > k x and k x > k y respectively. The region kx , ky <C 
where the results are inapplicable is not shown. The param- 
eter d is (a) d = 0.04, and (b) d = 0.13. Dashed line shows 
the mean field Fermi surface as given by Eq. Q. 



where n = (0,0), r 2 = (x 2 ,yi), Pi = (0, yi), and p 2 
(£2,2/2) (see Fig. |3]). Equation ([21~j) gives 



C (4) 



,4<Z 



2/1 (2/2 



Vi 



,_>(/ 



((2/2-2/i) 2 +£|) (2/i+£|) C 
(x 2 + y 2 2 ) d xl d 



(A2) 



We expect the main contribution to come from the re- 
gion 2/1,2/2^2/1 ^ x 2 it is convenient to introduce new 
variables, y 2 = £,x 2 , 2/1 = v£ x 2 ■ To isolate the leading sin- 
gularity in Eq. (IA1|) we expand the cumulant in Eq. (|A2[) 
for £ < 1 



C< 4 > « -2da 4d C 2 ^ 4 ^ 1 - 2<i (l - V ) 
This expansion holds for any d < 1/2. 
£< 4 > = -2di( Ja d ) 4 T(3 - 4d) 



1 — 2d 



(A3) 



/>00 />1 

/ d£ / ^77- 











3-4d 



(A4) 



The integral in Eq. (|A4|) is not convergent at the upper 
limit. This is an artifact of the approximation (|A3|) which 
is not justified for £ > 1. To overcome this we differenti- 
ate equation (|A1|) twice with respect to the parameter cv 
defined by Eq. (fHj) . 



9a 2 



2diT(5 - 4d)(Ja d ) 4 (-i(w - fc a )) 4 



/>00 r-\ 

< I dH drj 
'0 Jo 



£3-4d^l-2d^ _ ? ^l-2d 



(£ + a) 



5-4rf 



(A5) 



The integrations in Eq. (|A5[) are easily done with the 
result 



9 2 E (4) 



2dia" 1 r 2 (2 - 2d)- 



(Ja 



d\4 



(A6) 



9a 2 v ' (-i(uj - k y )) 3 - 4d 

Integrating Eq. (|A6|) back we finally get 

£ (4) = - 2dir 2 (2 - 2d)(Ja d ) 4 {~i{Lu - k y )) M - 3 

x (a log a + Aa + B) (A7) 

with A, -B constants. In the last equation the most sin- 
gular term is presented in Eq. (flT 



Fourth order contributions 



In the fourth order the general expression ()16|) takes 
the form 



£(?)(&) =iJ 4 / dx 2 e^ w - k ^ X2 dy 2 e 1(u - k ^ V2 
Jo Jo 

/ dy 1 C^(ri,r 2 ; Pl ,p 2 ), (Al) 
Jo 



2. Order J b 

In this case the general expression (fT6| with n = 3 
reduces to 

/>oo />oo 

E£ 6) (fc) =i(-i Ja d f dx 3 e [ ^- k ^ X3 dy 3 e 1( - u - k ^ V3 
Jo Jo 

r%3 fV3 rV2 

x dx 2 dy 2 / dyi 
Jo Jo Jo 

ABC-B-C + l (Ag) 



2/i (2/2 - 2/i) 2d (2/3 ~ 2/2) 



\2d 
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where 



self energy correction, 



A = 



(4 + y 2 2 ) d \xl + (y 3 - yi )r 



(A9) 



2 If/ 



B = 



and 



C 



[(x 3 - Xl ) 2 + (y 2 - yi) 
(x 3 - Xi) 2d 

x Kga - ^i) 2 + (y 3 - 2/ 2 ) 2 ] d 
X [(2/3 - 2/1) 2 + (a* - zi) 2 ] d 

\x\ +yl\ d \x\ + (y 2 - yi f\d 



(A10) 



(All) 



\ Xl ny 2 +xl\* 
We write 

ABC-5-C+1 = (A- 1)BC+(B-1)(C-1) . (A12) 



It is apparent form Eq. (|A12[) that the leading singularity 
originates from the first term, (A — 1)BC « A — 1 ~ 
— 2d(yi — 0)(j/3 —2/2)- Introducing new variables as in the 
Sec. IIHI and integrating over Xi s we obtain 



E (6) = 2di 



. ( Ja d fT{b - M) 



el— 2d t— 2d el 
SI S2 S3 



(a + 6 + 6+6) 



5-6ci 



(A13) 



Here again, the integrals are divergent on the upper limit. 
Similarly to the previous section we differentiate it once 
with respect to the variable a introduced in Eq. (|14[) in 
order to isolate the leading logarithmic singularity, 



da 



x (Jammed) 



(-i( U -ky))*-°* 

t\-2d c -2d c \ 
SI s"2 S3 



1 n % + \4 + v^ - <a") 



The remaining integrals are easily evaluated. The sub- 
sequent integration over a restores the singularity in the 



S (6) ^ d -KJa d fr(l~2d)T 2 (2-2d) 

(-i( W - fc,))5-6^ 

x(loga + C), (A15) 

where C is an integration constant. The singular part of 
Eq. (|A15[) is given by Eq. 



Appendix B: Leading singularities at the shadow 
mass shell, cj ~ fc H to the fourth order. 

In this appendix we evaluate the singular contributions 
to the Green function at the shadow mass shell, u ~ k y 
to fourth order in the coupling constant. We start with 
the expression (|A1[) introduced in App. IA 11 In contrast 
to the discussion in App. IA II we anticipate the singular- 
ity at u — k y to come from the region y<x 3> x 2 , and 
introduce new variables accordingly, x 2 = 6/21 J/i = VV2- 
Performing integration over y 2 we obtain 



EW =\{Ja d f 



m\ r(3--W),-«(i-,)-« 



Jo [(-i((w - fc„) + (w - fc x )£) 

(1 - v ) 2 + e\ d \v 2 + e\ d 



\i + e\ d \Z\ 2d 



1 



3-4d 



(Bl) 



We notice that the singularity at u) ~ fc y comes from the 
region of small £. Therefore we keep only the first term 
in the square brackets in Eq. (|B1[) . After performing 
remaining integrations over ^s we obtain 



(4) i( Ja d ) 4 r(l - 2rf)r(2 - 2rf)a 
(-i(w - fc.)) 3 - 4d 



2-2d 



(B2) 



We stress that contrary to the mass shell singularities 
discussed in App. [A] where it was important to compute 
the self-energy, at the shadow mass shell it is enough to 
consider the Green function itself. 
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